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Abstract. We explore the formation and evolution of the black hole X-ray 
binary system M33 X-7. In particular, we examine whether accounting for sys- 
tematic errors in the stellar parameters inherent to single star models, as well as 
the uncertainty in the distance to M33, can explain the discrepancy between the 
observed and expected luminosity of the 70 companion star. Our analysis 
assumes no prior interactions between the companion star and the black hole 
progenitor. We use four different stellar evolution codes, modified to include a 
variety of current stellar wind prescriptions. For the models satisfying the obser- 
vational constraints on the donor star's effective temperature and luminosity, we 
recalculate the black hole mass, the orbital separation, and the mean X-ray lu- 
minosity. Our best model, satisfying simultaneously all observational constraints 
except the observationally inferred companion mass, consists of a 13 black 
hole and a ~ 54 companion star. We conclude that a star with the observed 
mass and luminosity can not be explained via single star evolution models, and 
that a prior interaction between the companion star and the black hole progen- 
itor should be taken into account. 



1. Introduction 

M33 X-7 is the first stellar- mass black h o le (B H) discovered in an eclipsing 
X-ray binary system (XRB). lOrosz et al. ( 20071 ) constrained the mass to be 



15.65 Mq for the BH and ~ 70 Mq for t he 0-star cornpanio n. The orbital 
period {P) of ~ 3.45 days was determined bv lPietsch et al. ( 20061 ). The distance 



has been measured to be between ~ 840 kpc (jOrosz et al 



20071 ) and ~ 960 kpc 



(jBonanos et al.l l2006l ). Table [D and E] summarize other system's parameters 



relevant to our analysis. 

Table 1. Observed parameters for M33 X-7. 



Mdon(M0) 


70.0 ± 6.9 


Mbh(Mq) 


15.65 ± 1.45 


-RdonCR©) 


19.6 ± 0.9 


a(R0) 


42.4 ± 1.5 


P(days) 


3.45301 ± 0.00002 




74.6 ± 1.0 


/(Mbh)(M0) 


0.46 ± 0.08 




5% to 40% 


Spectral Type 


07 III to 08 III 


roff(K) 


35000 ± 1000 



Orosz and collaborators, private communication 2008 



It is challenging to make such a tight system with such a massive BH 
through standard BH-XRB formation channels. The current idea is that the BH 
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Table 2. Luminosity a nd X -ray luminosit y of M 33 X-7 for the distances 
adopted by lOrosz et all ()2007[ ) and Bonanos et all (120 061. The mean value 
for t he X-ray luminos ity according to .Orosz et al.l (|2007i) has been calculated 
from lLiu et all (|2008f) . 



Orosz et al. (2007) 

d(kpc) 840 ± 20 

log(Ldon/L0) 5.72 ± 0.07 

Lx(erg/s) (2.01 ± 0.48) xlO^^ 



Bonanos et al. (2006) 

(i(kpc) 964 ± 54 

log(Ldo„/L0) 5.84 ± 0.09 

Lx(erg/s) (2.65 ± 0.70) xlO^^ 



currently accretes mass from the wind of the 0-star companion. Mass transfer 
through Roche-lobe overflow (RLO) is excluded by the extreme mass ratio {q 
~4.5) since RLO would be dynami cally unstable and lead to a rapid merger 
of the two stars. Orosz et al.l ( 20071 ) also reported that the companion star is 
less luminous than what is expected from single star models. Because of the in- 
stability of mass transfer, this discrepancy cannot be explained by a past RLO 
phase from the 0-star to the BH (mass transfer on the thermal time scale of the 
mass donor in particular is known to produce underluminous donor stars). 



2. Modelling the System 

The goal of our analysis is to explore whether the observed system parameters, 
in particular the underluminous nature of the BH companion, can be explained 
via single star models when systematic errors in the star's parameters inherent to 
stellar evolution codes, as well as uncertainties in the distance to M33, are taken 
into account. We evolve a range of single star models accounting for mass loss 
through stellar winds, and extract the mass (Mdon) and radius (i?don) when the 
model matches the observed effective tempe rature (Tpff) and lu minosity (-Ldon)- 
We consider both the luminosity reported bv lOrosz et al. ( 20071 ) for a distance to 



M33 of ~ 84 kpc and a lurn inosity rescaled to the distance to M33 of ~ 964 kpc 
reported bv iBonanos et al.l (|2006) (Table [2]). Given the mass of the secondary 
when the model matches the observed Teg and L^^n we recalculate the mass of 
the BH (Mbh) from the observed mass function (/(Mbh)) and inclination [i). 
Next we calculate the orbital separation (a), the secondary's Roche-lobe radius, 
and the X-ray luminosity (i^x)- We restrict our investigation to those stars that 
are currently not filling their Roche-lobe. 

2.1. Uncertainties in Modeling Massive Stars 

The main uncertainties in the single star evolution models are the wind mass- 
loss rates and systematic errors in radius, luminosity and effective temperature 
inherent to stellar evolution codes (see Figure [1] for an example). We therefore 
use four different c odes modified to include an up-to- date spectrum of stellar 
wind prescriptions ( Hurlev et al.[2000l : Vink et 311120011 . and references therein). 



We then take into account the difference in the calculated radii, effective tem- 
peratures and luminosities for every considered single star model. To account 
for the uncertainty in the observed metallicity (Z), we consider models with 
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Figure 1 . 
in units of 
lution codes: 



Radial evolution of a 50 star as a function of time 
the main-seq uence li f e tim e for four different stellar evo- 
200l . 



_SSE (Hur ley et al 



, , _ EZ (iPaxtonI l200l. Eggleton 

jEggleton fc Kisele va-Eggleton 1 12002 ^ and EV ('Pod siadlowski et al. I2002bt 
llvano va et al. 200j|). At the end of the Main Sequence the systematic uncer- 
tainty in the radius can reach 5 . 



metallicities ranging from Z = 0.05 Zq to Z 
we present results only for Z = 0.2 Zq. 



0.5 Zq, but in the present paper 



2.2. The Mass of the O-Star From Single Stars Models at Z = 0.004 

The current observational constraints on the companion's radius, luminosity 
and effective temperature, narrow the range of possible masses that can be 
considered. Figure [2] shows HR diagram for different initial secondary masses 
calculated with the EZ and SSE stellar evolution codes. The SSE tracks cross 
the observed range of T^s one time. The EZ tracks cross the observed Tes three 
times, with the third crossing being in the Hertzsprung gap (HG). Given the low 
probability of observing a star in the HG we do not consider the third crossing 
any further. Using a distance to M33 of ~ 840 kpc, the observed luminosity 
and effective temperature constraints can be satisfied by masses up to 49 Mq for 
the l*** crossing found with both SSE and EZ, while the constraints can only be 
satisfied by masses up to 46 Mq for the 2"'^ crossing found with EZ. Assuming a 
distance of ~ 960 kpc , the constraints can be satisfied by masses up to 59 Mq 
for the I'** crossing found with both SSE and EZ, while the constraints can only 
be satisfied by masses up to 56 Mq for the 2"^^ crossing found with EZ. 



Results and Conclusions 



The main results of our analysis are summarized in Figure [3] and in Table [3?] 
If we assume a distance to M33 of ~ 840 kpc, the maximum value for the X- 
ray luminosity calculated with EZ is ~ 1.2 x 10^*^ erg/s, which is in agreement 
with the observed value at this distance to within a factor of ~ 1.2. On the 
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Figure 2. HR diagram for different initial ZAMS secondary masses and 
metallicity Z = 0.004 calculated with the EZ and SSE stellar evolution codes. 
The grey area represents the observational constraints on Toff, while the 
hatched areas represent the observational constraints on Ldon for distances 
of 840 and 960 kpc. 



other hand, adopting a distance of ~ 960 kpc, the maximum X-ray luminosity 
calculated with EZ is ~ 2.65 x lO^^erg/s which is in agreement within the 
error bars with the observed X-ray luminosity at this distance. As far as the 
0-star mass is concerned, the highest initial companion masses satisfying the 
effective temperature and luminosity constraints at ~ 840 kpc are 50 and 
46 M0 for SSE and EZ respectively. On the other hand, assuming a distance 
of ~ 960 kpc, we can consider masses up to 59 M© and 56 M© for SSE and EZ, 
respectively. Therefore, a zero-age main-sequence star (ZAMS) of 56 M© is the 
most massive single star model that satisfies simultaneously the observational 
constraints on the effective temperature, donor luminosity and X-ray luminosity. 
The calculated masses for this model for the 0-star and the BH are ~ 54 M© and 
~ 13 M©, respectively. We conclude that a companion star of 70 M© with the 
observed luminosity (log(L(ion/L©)~ 5.72) can not be explained via single star 
evolution models. Accounting for systematic errors in radius, luminosity and 
effective temperature inherent to stellar evolution codes, as well as uncertainties 
in the distance to M33 does not solve this discrepancy. 



4. Present and Future Work 

It is important to emphasize that our analysis treats the companion star as a 
single star, and does not account for any interaction between the companion star 
and the BH progenitor. Since single star models can not explain the observed 
properties of the 0-star, the system likely underwent a mass transfer phase 
during whi ch the BH progeii i tor fil led its Roche-lobe and transferred mass to 
the 0-star. Braun &: Langer ( 19951 ) shown that, depending on the efficiency of 
semi-convective mixing, such a mass transfer phase can lead to underluminous 
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Figure 3. X-ray luminosity as a function of present day donor mass for Z 
= 0.004, according to single stars evolution models, for distances of ~ 840 
(panel a) and ^ 960 kpc (panel b) . Accounting for the derived upper limits on 
the initial secondary mass for the 2"'' crossing of the observed Tcff constraints 
(see Figure [2]), we could only consider ZAMS masses up to 56 M©. The dashed 
horizontal line represents the observed X-ray luminosity. 



Table 3. Donor mass, BH mass, orbital separation, and maximum X-ray 
luminosity calculated for the highest initial donor mass satisfying the Toff and 
Ldon constraints (50 M© and 46 M© for SSE and EZ respectively for a distance 
of ^ 840 kpc, 59 Mq and 56 M© for SSE and EZ respectively for a distance 
of - 960 kpc). 



Code SSE EZ SSE EZ 

rf(kpc)-Crossing 840 - l*** 840 - 2"'' 964 - 1"* 964 - 2'"^ 

Mdon initial (M©) 50 46 59 56 

Mdon @ present (M©) 47.45 ± 0.03 44.71 ± 0.01 54.9 ± 0.2 53.99 ± 0.01 

Mbh @ present (Mq) 12.2 ± 0.8 11.8 ± 0.8 13.4 ± 0.9 13.2 ± 0.2 

a @ present (R©) 37.5 ± 0.2 36.9 ± 0.2 39.3 ± 0.2 39.1 ± 0.2 

Max Lx (erg/s) ~ 1.66 x lO^^ - 1.18 x lO^^ - 5.39 x lO^^ - 2.65 x lO^^ 



secondary stars (the so-called non-rejuvenation scenario). Therefore, we are cur- 
rently studying mass transfer sequences to explore this possibility. This analysis 
will also yield theoretical constraints on the age of the system, the binary com- 
ponent masses, the orbital separation right before the supernova explosion that 
formed the BH, and the mass lost and any possible natal kick imparted to the 
BH during the supernova explosion. 
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